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Abstract. The contribution of the pair annihilation pro- 
cess in relativistic electron-positron jets to the 7— ray 
emission of blazars is calculated. Under the same as- 

' sumptions as for the calculation of the yield of the in- 

. verse Compton scattered photons by relativistic particles 
in the jet (Dermer and Schlickeiser 1993) we calculate the 
emerging pair annihilation radiation taking into account 

■ all spectral broadening effects due to the energy spectra of 
the annihilating particles and the bulk motion of the jet. It 
is shown that in spectral appearance the time-integrated 
pair annihilation spectrum appears almost like the well- 
known 7-ray spectrum from decaying 7r°-mesons at rest, 
since it yields a broad bumpy feature located between 50 
and 100 MeV. We also demonstrate that for pair densi- 
ties greater 10 9 cm~ 3 in the jet the annihilation radiation 
will dominate the inverse Compton radiation, and due to 

, its bump-like spectrum indeed may explain reported spec- 
tral bumps at MeV energies. The refined treatment of the 
inverse Compton radiation including Klcin-Nishina cross 
section corrections, the low-energy cutoff in the radiating 
relativistic particle distribution function and a realistic 
greybody accretion disk target photon spectrum leads to 
spectral breaks of the inverse Compton emission in the 
MeV energy range with a change in spectral index Aa 
larger than 0.5 as detected in PKS 0528+134 and 3C273. 

Key words: Gamma rays - galaxies: nuclei - galaxies: 
active — black hole physics 



1. Introduction 

Until now 45 extragalactic 7-ray sources (status October 
1994) have been clearly detected by the EGRET exper- 
iment aboard the Compton Observatory (Kanbach 1994, 
Thompson 1994). Several of them are also seen by the 
COMPTEL and OSSE instruments (Kurfess 1994). The 
broadband high energy spectra extend over four decades 
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in energy and exhibit quite pronounced spectral forms and 
breaks. Apart from the Large Magellanic Cloud all extra- 
galactic 7-ray sources are radio-loud, have a compact core 
and are known as highly variable sources. A large fraction 
of these sources exhibit apparent superluminal motion. 

Because of its rapid variability, its high compactness, 
lack of radio-quiet sources in the sample and high pres- 
ence of superluminal motion sources it is generally agreed 
(Blandford 1993, Dermer 1993, Dermer and Schlickeiser 
1992, Dermer & Gehrels 1995) that the 7-ray emission 
originates in strongly beamed sources, in accord with the 
relativistic beaming hypothesis that has served well as the 
baseline model for the central engines of active galactic nu- 
clei (Scheuer and Readhead 1979, Blandford and Kdnigl 
1979). Dermer and Schlickeiser (1993) - hereinafter re- 
ferred to as DS - have demonstrated that at hard X-ray 
and 7-ray frequencies inverse Compton scattering of accre- 
tion disk photons by relativistic electrons and positrons in 
an outflow ing plasma with relativistic bulk motion (bulk 
Lorentz factor Y ~ 5 — 10) dominates over other radiation 
mechanisms. This 7-ray production process succesfully ex- 
plains the key observational features as (i) the peak of the 
bolomctric luminosity at 7-ray frequencies (ii) the 7-ray 
variability time scale < few days, (iii) the softening of the 
spectra between the hard X-ray and medium 7-ray regime, 
(iv) the high percentage of superluminal sources in the 
sample, (v) the generation of TeV emission for favourably 
aligned observers. 

It has been suggested (Henri, Pelletier & Roland 1993) 
that pair annihilation in relativistic electron-positron jets 
could be of considerable importance to such jet models. 
The purpose of this study is thus to extend the DS-model 
by including the 7-ray production from pair annihilation 
radiation from relativistic jets. Under the same assump- 
tions as for the calculation of the inverse Compton scat- 
tering from the jet we will determine the yield of anni- 
hilation radiation. Of special interest will be the mod- 
elling of spectral breaks with a change in spectral index 
Aa larger than 0.5 as detected in PKS 0528+134, Aa = 
1.16(+0.23,-0.11) and 3C273, Aa = 0.66(+0.13, -0.15). 
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In the original DS approach spectral breaks up to Aa < 
0.5 were obtained by incomplete inverse Compton cooling 
of an injected power law during the time evolution of the 
7-ray flare. Taking into account all spectral broadening 
effects due to the energy spectra of the annihilating parti- 
cles and the bulk motion of the jet, we show here that the 
contribution from this radiation component peaks near a 
few McV, and if it dominates the inverse Compton com- 
ponent it may provide an explanation for the reported 
spectral maxima in this energy range. 



many superluminal motion sources in the sample of 7-ray 
blazars. At these heights the pair annihilation radiation 
freely escapes from the system and does not initiate a 
pair cascade (Lightman and Zdziarski 1987, Zdziarski et 
al. 1990, Protheroe et al. 1992) with plenty production 
of 511 keV annihilation line radiation from the cascade- 
cooled cold electron-positron pairs that is not observed 
from these objects (Kurfess 1994). As in DS we start wit 
h initial relativistic power law energy spectra for the ra- 
diating particles 



2. Jet model for the 7-ray emission 

The jet model for the 7-ray emission in AGNs of DS is il- 
lustrated in Figure 1. Relativistic electrons and positrons, 
with isotropic energy distribution functions n±('j; z) in 
the comoving frame at height z above the accretion disk, 
f low outward along the disk symmetry axis with bulk 
Lorentz factor T and bulk velocity /?rc from an initial 
height Zi. Photons produced by the accretion disk pass 
through the jet and are Compton upscattered to 7-ray 
energi es km e c 2 (see DS for details). The relativistic elec- 
trons and positrons fill not the whole jet but only an (as- 
sumed) spherical blob of radius Rb- We follow the nota- 
tion of DS: energies and angles without asterisks denote 
values in the comoving fluid frame, while these quantities 
with asterisks denote values in the stationary frame K of 
the accretion disk. We do not include cosmological effects 
here. 



e -e -Plasma 




Accretion disk 
Fig. 1: Jet model for the 7-ray emission in AGNs 



«±(7 



= N ±1± S l«7i± <7± <72±, (1) 



and follow the time evolution of the energy spectra 
under the various energy loss processes. These spectra are 
used to calculate the yield of annihilation radiation at var- 
ious heights. We start our analysis with the calculation of 
the annihilation radiation emissivity. 

3. Pair annihilation radiation from the electron- 
positron blob 

The general calculation of the yield of pair annihilation 
radiation from relativistic plasmas has been treated by 
Svensson (1982). In the case of an ultrarelativistic pair 
plasma, characterised by Eq. (1), this treatment can be 
considerably simplified by using the differential pair anni- 
hilation cross section 



da 



rel 



Y 2 
I cm 



4x 

which is valid for -f c 



1 



3> 1. The subscript cm denotes 
quantities measured in the center-of-momentum frame 



and x 



ing to Svensson's formulae (55) 



. Then, following the calculation lcad- 
(57) yields 



da 

dk, 



1 



Trr„ 



rel 



+ 



lira V(l+ +7-) 2 -47, 

2 2 

7cm \flcra ~i~ c +7cm 7cm \/7cm 

where 



2 

cm 



+ 



(3) 



To avoid heavy attenuation of the produced GeV 7- 
rays by pair production interactions with the accretion 
disk photons Dermer and Schlickeiser (1994) have pointed 
out that the initial height of the blob has to be larger 
than z p = lO 15 (L46/05o) 1 ^ 2 cm where L^q denotes the ac- 
cretion disk luminosity in units of 10 46 erg/s and 650 the 
accretion disk temperature in units of 50 eV. Such a high 
location of the blob is also consistent with the presence of 



1 cm 



(k- lT f-l, (4) 

fc(7++7--*0, (5) 
k is the photon energy normalized to the electron rest 
energy, the subscript ± denotes quantities rcfering to the 
motion of the positron/electron as measured in the la- 
bor frame (the rest frame of the plasma blob) and 7* m is 
the maximum center-of-momentum 7-factor in the case of 
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which a photon of energy k can be produced, which is the 
notation used by Svensson (1982) — and 



dk Jrel 



7 2 7 2 



(v / (7++7-) 2 -4 7c 2 m + H + + H-) 



i u cm 



(6) 



(12) simply becomes k < 7+ + 7- which is equivalent 
to the condition 7* m > 0. 

Unfortunately, equation (6) is still not analytically in- 
tegrable, but for 71+, 71 _ £ 5, it may be very well used 
for numerical calculations of the integral 



where the bar denotes the averaging over collision an- 
gles between (3 + and (3 _ and 

7 cm ITm, 



7 *2 

I cm 



arcsin 



l c ± 



3/2 



— arsinhl ■ 



1 cm 

/c± 7c 



V la 



H± = 



111 

3 7 c * 



3 

cm 



C±7cm + 

for c± < 
for c± > 

for c± = 



(7) 



(8) 



Also the integration limits 7cm, 7^ which are deter- 
mined by kinematical restrictions may be simplified. As 
given by Svensson (1982) these integration limits are found 
to be 



~L _ min 
I cm 1 cm 



-(l + 7+7-[l-^-D 



r max * l 
I I cm 1 1 cm J 



where 



max 
I cm 



(1 + 7+7- [1 + /3+/3-]) 



(9a) 
(96) 

(10) 



is the center-of-momcntum 7-factor in the case of 
head-on collision which, for /3+, /?_ as 1, simplifies 
to 7c„ :l: » -\/7+7-- 7^™ is the minimum center-of- 
momentum Lorentz factor in the case of an overtaking 
collision. Writing 7* m = 7+7- — (k — j-)(k — 7+), we 
immediately see that 



1 U = 

I cm 

7™ x for (k > 7+ A k < 7_) or (fc < 7+ A k > 7_) 
7* m for (fc < 7+ A fc < 7_) or (fc > 7+ A fc > 7_) 

(11) 

Furthermore, we have to satisfy the condition 



72+ 

- / * 

71+ 



72+ 72- 

/ dl+ / 



07. 



da , 



f+ (7+) /- (7- ) v-^ (k, 7+ , 7_ 



(13) 



71- 



(Svensson 1982) where the brackets denote the av- 
eraging over the distribution functions of electrons and 
positrons. 

For further simplifications, we use the fact that in 
Eq. (6) only the upper limit is important and also 
the first term is negligible as compared to the values of 
H±{k,'~fc m ,7c m )- As mentioned by Svensson, for highly 
relativistic pair p lasmas, the angle averaged emissivity is 
sharply peaked at 7+ = k and 7_ = k where the pho- 
ton takes over the energy and momentum of the incoming 
positron or electron, respectively. We use this fact to ap- 
proximate v^- by constructing a function that satisfies 
the condition of these peak values and exhibits the same 
asymptotical behavior as found for vj^. 

The well-known <5-function approximation does not 
work well in the case of an energy distribution as given by 
eq. (1), particularly, it does not describe the lower energy 
tail of the annihilation spectrum (k < min {71+, 71-}) at 
all. The approximation scheme we derive here is much sim- 
pler than that given by Coppi and Blandford (1990) who 
still need numerical integrations to determine the disper- 
sion of u4| about k but whose results are — in contrast 
to ours — also valid for midly relativistic pair plasmas. 

The peak values of H± as given by eq. (6) are easily 
found by setting c± = -1 and 7^ = 7 * m = ^/t+tT: 



H±(k,jL) 



2 V 1 cm J 



(14) 



Considering the asymptotic behavior of v^, we see 
that for |fc— 7± | » 1 the photon is not emitted in the direc- 
tion of the incoming electron/positron, and so x does not 
approach 1 where the differential cross section (2) peaks. 
For this reason, we may set f3 cm = 1 in cq. (2) and in the 
following steps of calculation, which yields 



= 5 (7 + [l 



k < ku 
f-/?+]+7-[l 



(12) 



where ku is the maximum photon energy available by 
the annihilation of an e+-e~-pair characterized by 7+, 7_ 
(Svensson 1982). In the ultrarelativistic case, equation 



da 

^(fc,7+,7-) 



C±»l 



I ^ — T=F I 



(15) 



A quite simple function combining the characteristics 
described by eqs. (14) and (15) and therefore approximat- 



ing v% is 
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(7 
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cm ) 



1^-7+1 + ^ I* -7-1 + 
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Fig. 2: Asymptotic approximation (16) (transparent area) 
in comparison to the relativistical value of u4r 
calculated by eq. (6) (gray area) for k = 4 





Fig. 4: Same as Figs. 2 and 3 for k = 100 



distr ibution functions. Due to different determinations of 
u4? in the four integration areas defined by Eq. (11), the 
results look somewhat involved, but they are — after some 
manipulations — numerically quite simply evaluable. 

Expanding (16) in a series and performing the integral 
(13) — where we neglect the integration area given by 
the edge points (7+, 7+) = ( k , k ~ 7i±), (k, 7i±), (7i±, k - 
7i+),(7i+, fc -7i±),(7i+,*0 for integration over i +l l_ 7± 
— yields 



4> « 



e(k - 7i-)e(fc - 71+) 



7i_ - (min{fc,7 2 _})~ 



[5'i(min{fc,72+}, s + 1) - 5i(max{7i+, k - 71-}, s + 1)] 

7i- - (min{A:,72-}) 
s + 1 

• [5i(min{fc,7 2+ },s) - S'i(max{7i + , A: - 71-}, s) 
- /s{5i(min{fc,7 2+ }, s + 1) 



-Si(max{7i + , k - 7 X _}, s + 1)}] 



(17a) 



Fig. 3: Same as Fig. 2, but for k = 15 



We call it the asymptotic approximation. Figures 2 
- 4 show that Eq. (16) is a very useful approximation 
to the angle averaged emissivity from highly relativistic 
pair plasmas for k ^, min {71+, 71-}, but which is — in 
contrast to the (^-function approximation — still sensible 
for k < min{7 1+ ,7!_}. 

By means of a series expansion of the asymptotic ap- 
proximation, the integrations in eq. (13) are elementary 
for power-law distributions (1) and are also analytically 
executable for other physically relevant ultrarelativistic 



+ 6(72- - k)Q(k - 7l+ ) 



(max{fc,7i_|) s -7 2 - 
k s 



(Si(min{fc,7 2+ },s) - Si(~f 1+ ,s)) (17b) 



+ Q(k — 7i_)6(72+ - k) 



(min{fc,7 2 -}) s - 7 X _ S 



(5*2(72+, s + 1) - S* 2 (max{fc,7 1+ },s+ 1)) (17c) 



+ 6(72- - fc)6(72+ - k) 



7 2 - ^ (max{/c,7i-})~ 



■ k (52(72+, s + 2) - S*2(max{fc, 71+}, s + 2)) 
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7 -j s+1) -(max{fc ;7l ,})" ( - +1) 
+ IT! h 



■ (52(72+, s + 1) - S , 2 (max{fc,7i + }, s + 1) 
-k[S 2 {l2+,s + 2) - S 2 (max{fc, 7l+ }, a + 2)])^ (17d) 
+ [+-"] 



where [+ <-> — ] denotes the preceding terms (17 a) - 
(17 d) under transposition of the subscripts + <-> — , is 
the Heaviside function and 



S 1 ( 7 ,s)= 7 - s V +^fc- s ln( 7 ) (18a) 



n = 



S 2 (7,s)=7- S E 



1 



n=0 



n + s \ 7 



10 
10 2 
10 2 
10 2 
10 2 

„ 10 " 2 
™ M 2 

£ 10 2 

.O. 

10 

T3 

I 10 2 

V 

10 2 
10 2 
10 3 
10' 3 
10- 3 
10 3 




10000 



(186) 



Fig. 5: Pair annihilation spectra calculated by the asymp- 
totic approximation (17) (dashed) compared to numerical inte- 
grations (solid and dotted, respectively); parameters: 71+ = 8, 
71- = 15, 72- = 72+ = 10 4 



with the symbol 5 s ,n defined by 

^,jv = . , , _ r (19) 
I for s f IN 

and k± = k ± ^ . Of course, we must not replace n± by 
fc in eqs. (18), because the series Si and S2 would diverge 
then. The series Si is easily evaluable by summing sev- 
eral (> s) terms and replacing the rem aining part of the 
sum by the integral whose upper Riemann sum it is, while 
S2 may be evaluated by subtracting a geometrical series 
whose value is known, with the remaining series converg- 
ing very quickly. 

Doing so, we see that Eq. (17) provides a very accurate 
approximation to the pair annihilation spectrum. In Fig. 
5, this approximation is compared to the results of nu- 
merical integrations over the exact solution for vjfi^ given 
by Svensson (1982). Especially the peak structure and the 
high energy tail of the spectrum are very well described 
by the asymptotic approximation, while also the low en- 
ergy tail is sensibly approximated, considering that this 
part of the spectrum — as compared to spectra produ ced 
by other radiation mechanisms resulting mainly in power 
laws — will not be very important. Note that the different 
absolute values of the peak emissivity for different spectral 
indices is caused by the fact that for simplicity we used 
not normalized distribution functions with N± = 1. 



4. Particle spectra evolution 

We now consider the evolution of an ensemble of particles 
that is injected at height z% with a power law distribu- 
tion in energy given by Eq. (1). We first treat the effects 
of inverse-Compton, synchrotron bremsstrahlung and pair 
annihilation losses seperately in order to achieve an esti- 
mate for typical timescales for the deviation from an ini- 
tial power law distribution caused by each one of these 
processes. 

The typical peak structure of annihilation radiation is 
mainly caused by the lower energy end of particle distri- 
butions. Therefore we show that energy loss mechanisms 
most effectivly cooling high energy particles (7+ » 7 i±) 
may not seriously influence the shape of the annihilation 
spectrum during the evolution of the particle spectrum. 
So, the typical timescale for low energy particle spectra 
evolution may be significantly longer than the timescale 
for energy losses by a single particle. 

Note that this does not represent a realistic calculation 
of the particle spectra evolution — this will be treated in 
future work — , but gives only a rough estimate of the 
typical timescale on which the assumption of a power law 
distribution of particle energies is justified. 
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4-1. Inverse Compton losses 



{.2. Synchrotron losses 



The effect of inverse Compton losses on particle energy 
spectra in the jet has been treated in DS. There, a typical 
particle energy is defined by 



7o 



2488 £ 



/? r r 3 (l + /3 r ) 2 



(20) 



where the tilde denotes the normalization to gravita- 
tional radii — z = , R„ = — and £ a = It 

t^g a c ^edd 

is shown that particles with energy greater than are 
suffering heavy losses while the spectral shape of particles 
with energy j± < 700 remains nearly unaffected by inverse 
Compton losses. This spectral evolution which, of course, 
does not depend on the particles' charge is shown in Fig. 
6. 




10000 



520 1 



100000 1000000 



Fig. 6: Evolution of a power law particle spectrum, affected 
by inverse-Compton losses. Parameters: Zi = 500, 7i±i = 10, 
72±l = 10 6 , s = 2.5, M 8 = 1, £0 = 0.1. 



From Eq. (20) and Fig. 6, we conclude that the shape 
of the low energy particle spectrum will not be altered 
substantially by inverse Compton losses, and therefore the 
pair annihilation peak structure will not strongly depend 
on particle spectra evolution caused by the inverse Comp- 
ton process. 



In order to investigate the effect of synchrotron losses we 
use the single particle energy loss rate due to the syn- 
chrotron process 



-mc 2 7± 



4 £ 2 R2 2 



(21) 



- where B is the magnetic field strength which we 
assume to be approximately constant — combined with 
the blob's equation of motion 



z(t*) = Zi + 



c(3rt* 
R„ 



(22) 



and the time-dilation St* = TSt. This yields 



dz 



where 



Sy 



4 a T B 2 R g 2 

3 87TTOC 2 T% 7± 



-A"7±- (23) 



K = 5 • 10~ 7 



Mr 



The differential equation (23) is solved by 

l±\z) = l±l + K(z-Zi) 
which is valid in the bounds 

7i± := (li±i + K[z - Zi]) 1 < 7± 
7± < {I2L + K \z - 1 = : 72±- 



(24) 



(24) 



(25) 



Considering conservation of particle number and in- 
serting eq. (1) for the initial state leads to the resulting 
spectrum 



« ± (7±^)-^±7 ± 2 (7± 1 - K[z-Zi}Y 



(26) 



From Fig. 7 we see that the particle spectrum influ- 
enced by synchrotron losses is well represented by a power 
law with the original spectral index between the cutoffs 
given by Eq. (25). But in contrast to the case of inverse 
Compton losses also the lower energy cutoff may be con- 
siderably lowered and the higher energy cutoff does not 
approach a fixed limit if the magnetic field strength re- 
mains constant along the jet, as we assumed. 

A rough estimate for the blob's distance z from the 
accretion disk plane at which the lower energy cutoff 
is significantly influenced by synchrotron losses may be 
achieved by setting ~/i±(zo) = yielding 



z t + ^—~z t + 2-W-^—(^) . (27) 
1i±iK 7i ±i M 8 \GJ 
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4-3. Bremsstrahlung losses 

An estimation of particle energy losses due to relativis- 
tic e + -e~ and e~-e~ bremsstrahlung is obtained using 
the formalism similar to the one proposed by Dermer & 
Liang (1989) for power law particle spectra. In contrast to 
their calculation, we include the Doppler boosting of the 
bremsstrahlung photon assuming that it is emitted in the 
direction of one of the incoming particles. Furthermore, we 
use the ultrarelativistic limit of the bremsstrahlung cross 
section derived by Baier et al. (1967). 

Results of this treatment are shown in Fig. 8 which 
shows that the single particle energy loss rate is roughly 
proportional to the particle energy so that the energy loss 
timescale remains roughly constant. 

As we also see from our numerical calculations, there 
is no very significant dependence of the energy loss time 
scale on the lower cutt-off of the particle distributions. 
Considering the lower energy end of the particle spectra, 
we find a typical energy loss timescale of the order 



Tbr = 5 • 10 s ■ n 



- where nn = 
blob distance of 



mi 



(28) 

which corresponds to a 



4-4- Thermalization 

The thermalization timescale for a plasma described by 
a power law in particle energy is again estimated using 
the formalism of Dermer & Liang (1989) to calculate the 
single particle energy loss and gain, respectively, due to 
M0ller and Bhabha scattering. The numerical results for 
M0ller scattering for several values of the power law spec- 
tral index s are shown in Fig. 9; the difcrcncc between 
M0ller and Bhabha scattering turns out to be negligible. 
We only consider the timscales under the initial condi- 
tions; synchrotron, inverse Compton and bremsstrahlung 
losses, of course, will alter the upper cut-off of the particle 
distributions and hence, they will lengthen the thermal- 
ization timescale for particles with lower energy. 

Obviously, the lower energy end of the particle spec- 
trum is most seriously affected by energy gain due to 
M0ller scattering. From our numerical calculations wc find 
the timescale for this energy exchange process as 



T th = 20s 



7i- 7i+ 



nn 



(30) 



corresponding to a distance 



zbr 



h + 10 3 



Ti 



n n M 8 



(29). 



zf =~ Zl + 0.4 



£1 71 + 71- 

mi m s 



(31) 
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Fig 9: Timescales for single particle energy loss or gain, 
respectively due to M0ller scattering; n± — 
10 11 cm- 3 



4-5. Pair annihilation losses 

In the general case, the particle loss rate can be written 

as 



dt 

72 =p 



PA 



In order to estimate the timescale for the evolution 
of initial power law particle spectra we may replace 
n T (7 T ,z) by its initial value given by eq. (1). Then, the 
integral in eq. (33) is analytically solvable yielding 



dz 



PA 



-NHl ± ,~z)^^ (34) 
7it7± 



where, for spectral indices 2 < s < 3, P may be 
roughly estimated by 



P = 



(s - l)A7rr e 2 



1 



,r Pr V M47 ± 7t T ] + --2 

in-12 M $ 3 
10 ——cm 




dj T n T (7 T ,z) wcr(7 + ,7_) 



(32) 



7l T 



which, in the ultrarelativistic limit and for a homoge- 
neous blob, becomes 
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Fig. 10: Evolution of a power law particle energy spectrum 
influenced exclusively by pair annihilation losses 
in the case of an electron-dominated pair plasma. 
Parameters: n~ = 10 11 crrT A >• n+, z% = 500, 
7i±< = 10, 72±« = 10 6 , s = 2.5, M 8 = 1, F = 10. 



(35) 



dz 
72 T 



= -N ± ( 7± ,z)—-Z-- 
pa 1 Pr 7± 



We now consider two special cases in which the dif- 
ferential equation (34) is easily solvable. First, we as- 
sume the pair plasma to be dominated by electrons (an 
(33) e~-p + -plasma), thus n_ 3> n+. Further neglecting geo- 
metrical effects on the extent of the plasma blob, we set 
where N ± (j±,z) is the spectral particle number and n - = cons t-i an d obtain 
we used Svensson's formula (19 b) for the angle averaged n _ P 
reaction rate va. N + (j+,z) = N + (~/+, Zi) e~^-"'+ {z ~ Zi) (36) 



/ 



7l=F 



d lT n ^z) ( M 4 7+7 _]- 2 ). 

7t 
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which is shown in Fig. 10. From Eq. (36) we see the 
typical evolution time scale to be 



1 



zo 



T7i_ 7i + 



A M s n_ [cm J J 



(37) 



Secondly, we assume that electrons and positrons have 
the same particle density and distribution functions and 
again neglect geometrical effects eq. (34) becomes 



= -N ± \ 1± ,~z)- 

PA V b 7it7± 



( dN±{ 1± ,z) 
^ dS 

where V& is the blob's volume. This is solved by 

JV ± ( 7±) 5) 



co + 



Pi 
V b 7l±7± 



(38) 



(39) 



where the constant Co can be estimated by N ±^ ± ■ 
Thus we write 



7V ± (7±,^) - 



1 _l_ Z-Zj ' 



(40) 




Fig. 11: Evolution of a power law particle energy spec- 
trum influenced by pair annihilation losses in 
the case of electrons and positrons having the 
same density and distribution function. Param- 
eters: n± = 10 11 cm -3 , 2% — 500, 7i±j = 10, 
72±i = 10 6 , s = 2.5, M s = 1, T = 10. 



This behavior which is, of course only self-consistent as 
long as the particle spectra may be well fitted by a power 
law with spectral index <~ s is shown in Fig. 1 1 . We derive 
the typical evolution timescale to be 



zo(7±) : = 



7? ± 



Pn±(zi) 



10 1 



7i±r 



M 8 n±(zi)[cm- 3 Y 



(41) 



where we used — ^±i£i) — 



4.6. Comparison of times coles 

We now compare the typical timescales for particle spec- 
tra evolution by the energy loss and exchange processes 
considered above. 

First, we see that the typical energy loss timescale for 
the inverse- Compton process is much shorter than for pair 
annihilation as well as for all the other processes consid- 
ered above, but it only influences particles with energies 
above 700- Thus, it's negligible in calculations to the yield 
of pair annihilation radiation. 

For the other processes, we estimated that 



zo- 



10 2 7l±Fl 

1U Ms nil 
-^q2 71 + 71- 



11 (I)" 



10 

71 ± M « 

10 3 Vl — 

nil M s 

0.4 r ^i+-n- 



for PA, n + = ra_ 
for PA, n + <C n_ 

for synchrotron 
for bremsstrahlung 
for thermalization 



(42) 



Comparing these timescales we see that for lower cut- 
offs £ 3 the effect of bremsstrahlung losses are of minor 
importance to the particle spectra evolution. Particularly 
for the lower energy end of the particle spectrum, the 
bremsstrahlung losses are overestimated by the cross sec- 
tion approximation used in our calculation. So, this con- 
straint will be eased when treating the bremsstrahlung 
emission process in more detail. But because of being 
much more efficient than bremsstrahlung emission energy 
gain due to M0ller and Bhabha scattering will compensate 
the according energy losses. 

As can be seen in Fig. 9, thermalization will affect the 
lower energy end of the particle spectra very seriously. But 
since it is no energy loss process and thus — by its nature 
— does not lead to a cooling of the plasma the latter will 
remain ultrarelativistic. This energy exchange process, of 
course, will alter the shape of the particle spectrum at low 
energies, but the mean energy of the particles remains con- 
stant, and the higher energy end of the particle spectrum 
will not suffer heavy energy losses. Thus, for a detailed 
discussion of the particle spectra evolution thermalization 
effects must be taken into account, but they do not seri- 
ously alter the shape of the time integrated annihilation 
spectrum. 



10 



M. Bottcher, R. Schlickeiser: Pair annihilation radiation from relativistic jets in 7-ray blazars 



Furthermore, we can derive a critical magnetic field 
strength for the dominance of synchrotron losses over pair 
annihilation losses on particle spectra evolution. This is 
given by 



E>cr 

~G 



10 5 ^i 



for n + = ra_ 
for n + <C n_ 



(43) 



If B is significantly lower than B cr we may really ne- 
glect the effect of synchrotron cooling on the particle en- 
ergy spectrum. 

5. Comparison to the inverse-Compton spectrum 

We now compare the yield of pair annihilation photons to 
the 7-ray spectrum produced by inverse-Compton scatter- 
ing of accretion disk photons by the high energy particles 
of the jet. As we pointed out in section 4.6, this process is 
only efficient for high energy particles 7 > 700, but works 
over a much shorter timescale than all the other processes 
occurring in the pair plasma. According to the long in- 
tegration times of the observing instruments we only see 
time integrated spectra of AGN. 

Thus, it is of interest to compare the time integrated 
inverse-Compton spectrum to the pair annihilation spec- 
trum integrated over the typical timescale cither for pair 
annihilation effects on the lower energy particles or for 
any other mechanism reducing the particle density of pairs 
near the lower energy cutoff. 



10* 



Pair annihilation 




10 



10 



10" 10' :0 

k = hv/(m a c 2 ) 



10 J 



Fig. 12: Time-integrated inverse-Compton and pair annihila- 
tion spectra; parameters in the text. 



For the following calculations, we assume this 
timescale to be of the order of z ~ 10 3 . Each plasma 
blob is assumed to be spherical with a radius of 1 11 cm 
and injected at Z{ = 10. The accretion disk radiates with 
Lq = 10 46 ^p, and the mass of the central black hole is 



10 8 M Q . The accretion disk radiatiation spectrum is mod- 
elled by a greybody distribution with a radial dependence 
of the tempearure according to the accretion disk model 
of Payne and Eardley (1977). The lower energy cutoff of 
the distribution functions is 71-t = 3, the higher cutoff is 
72-t = 10 , and we assume that both electron and positron 
density are equal . The spectral index of the distribution 
functions is s = 2. The observing angle with respect to 
the jet axis is 9 = 2°, and the blob moves outwards with a 
Lorentz factor of T = 10. These parameters are convenient 
to explain the 7- ray spectra of some 7-ray blazars by the 
inverse-Compton process. 

The inverse-Compton radiation spectrum is calculated 
using the full Klein-Nishina cross section in the head- 
on approximation. The resulting photon number spectra 
from inverse-Compton scattering and pair annihilation are 
shown in Fig. 12 where we see that for particle densities 
n 10 9 cm~ 3 the contribution of pair annihilation should 
lead to an observable effect in the spectra. The fiat inverse 
Compton spectrum at small photon energies is particu- 
larly noteworthy. 

We attribute this mainly to the presence of a low- 
energy cutoff at E c — 7 C mc 2 in the radiating particle dis- 
tribution function. From the work of Schlickeiser (1982) it 
is well known that such a low-energy cutoff indeed leeds to 
a flat scattered photon spectrum below the photon energy 



4e 7 c 2 , 



(44) 



where eo denotes the energy of the assumed monoener- 
getic target photon distribution function. (Corresponding 
calculations for extended target photon energy distribu- 
tions by Schlickeiser [1980] showed similar flattenings.) 
According to Eqs. (14) and Fig. 1 of Schlickeiser (1982) 
the inverse Compton power of scattered photons varies 
proportional to 



x(Ej) cx E. 



1-3 
2 

'7 



for E. t > E» 



(45a) 



where s is the power-law index above E c , while at energies 
£ o < Ej <C E® one obtains 



x(Ej) oc 



E~ 



(45b) 



indicating that indeed spectral breaks in the inverse 
Compton scattered photon spectrum by Aa = 1 — = 
around the intrinsic energy E° are implied. Due to the 
blueshift by the relativistically moving blob this spectral 
break occurs at 



4eo7^ 



with the Doppler factor 



^ = L96 (To^y)(l) 2 ^ ^ 



d<jj?)= r(i-/3iV) 



(47) 
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k c lies in the MeV energy band which is in agreement with 
the inverse Compton curves in Fig. 12. With parameter 
values as chosen here, on top of this break the contribution 
from the pair annihilation radiation is located. 

5.1. Comparison with observations 

It can be seen from Fig. 12, that the pair annihilation emis- 
sion in such a highly relativistic plasma does not result in 
a 511-keV peak. The annihilation radiation is blucshiftcd 
both, by the relativistic motion of the pairs in the blob 
rest frame, and by the Doppler boosting due to the bulk 
motion. The pair annihilation peaks near the photon en- 
ergy 

k max =liD(p*) (48), 

with the Doppler factor (47). 

For the parameters chosen in Fig. 12 we obtain k max = 
89.2, so that the annihilation feature is at energies of 
<~ 50 - 100 MeV. In this case the annihilation bump ap- 
pears almost like the well-known (e. g. Cavallo & Gould 
1971, Stecker 1973) 7-ray bump from decaying 7r°'s near 
70 MeV if they do not move relativistically. The detec- 
tion of enhanced bumpy emission near 70 MeV from an 
AGN would therefore be no unique clue on the presence 
of hadrons in the source, but may point to blucshifted 
annihilation radiation. 

However, for a sample of AGN sources seen under dif- 
ferent viewing angles [i* there will be quite some variation 
in the values of the parameters 71 , defining the low-energy 
cutoff of the relativistic electron and positron energy spec- 
trum, and the relativistic bulk speed of the emitting blob, 
r, from source to source. While the energy exchange pro- 
cesses discussed in section 4, will keep the low-energy 
turnover of the distribution function at mildly relativis- 
tic values, 71 > 3, the wide dispersion in the Doppler 
factor (47) values will imply a corresponding wide distri- 
bution in the value of the annihilation peak energy (48). 
Although we have made no detailed effort so far to opti- 
mize parameters to explain individual source spectra, it 
is evident that in principle the inclusion of the annihila- 
tion radiation is capable to explain emission maxima and 
spectral breaks at MeV energies. The occurrence of spec- 
tral maxima between 1 and 10 MeV for the sources 3C273 
and 3C279 has been reported by recent Comptel observa- 
tions (Williams et al. 1995). An even more drastic example 
is the detection of the MeV-blazar GRO J05 16-609 near 
PKS 0506-612/0522-611 by Bloemen et al. (1995). Our 
spectral modelling confirms the suggestion of Bloemen et 
al. that the MeV emission in this object represents the 
broad blue-shifted annihilation line. 

A possible test for the hypothesis that the emission 
maxima in the observed 7-ray energy spectra are due to 
blucshiftcd annihilation radiation is offered by the pre- 
dicted relation between the maximum photon energy (48) 
and the total 7-ray luminosity. Making the strong assump- 
tion, that apart from the Doppler factor all other intrinsic 



physical parameters of the sources are equal, the total pair 
annihilation luminosity will scale with the Doppler factor 
as L p oc D 3 . A similar Doppler factor dependence holds 
for the inverse-Compton luminosity (DS). Combined with 
Eq. (48) this implies the relation 



Ll /3 



(49), 



i.e. the strong Doppler boosting in bright objects also 
blueshifts the annihilation peak to high photon energies. 



10 



10° 



10' 
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Fig. 13: Relation between MeV-7-ray peak energy and total 7- 
ray luminosity for 33 AGN detected by EGRET. 



We have used the available compilation of blazar en- 
ergy spectra (Montigny et al. 1995) to estimate both the 7- 
ray luminosity L p and the maximum photon energy (cor- 
recting for cosmological redshift), and the result is shown 
in Fig. 13 in comparison with the predicted correlation 
(49). The presently available data are not in contradic- 
tion with the hypothesis (49). A more stringent test would 
be possible if more AGNs with luminosities in the range 
0.0005 — 0.03L 48 are detected with well-resolved 7-ray en- 
ergy spectra. 

6. Summary and conclusions 

In this work we have calculated the contribution of the 
pair annihilation process in relativistic electron-positron 
jets to the 7— ray emission. Under the same assumptions 
as for the calculation of the yield of the inverse Compton 
scattered photons by relativistic particles in the jet (Der- 
mer and Schlickeiser 1993) we calculated the emerging 
pair annihilation radiation taking into account all spec- 
tral broadening effects due to the energy spectra of the 
annihilating particles and the bulk motion of the jet. We 
started from the pair annihilation reaction rate in rela- 
tivistic plasmas (Svenson 1982) and developed a new asym 
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ptotic approximation appropriate for relativistic pair en- 
ergies which provides a very accurate approximation of 
the pair annihilation spectrum for all photon energies. 

The evolution of the particle energy spectra with time 
has been taken into account consistingly, discussing in de- 
tail the role of inverse Compton losses, synchrotron losses, 
relativistic bremsstrahlung losses, pair annihilation losses 
and M0ller and Bhabha thermalization. For typical blazar 
parameters two particle-number conserving loss processes 
dominate at large (Lorentz factor 7 > 10 3 ) particle en- 
ergies (inverse Compton losses) and at very low (Lorentz 
factor 7 < 30) energies (thermalization), while the catas- 
tropic pair annihilation los ses reduce the number of radi- 
ating particles thus avoiding the " dead-electron-problem" 
(Begelman et al. 1984). We demonstrated that in spectral 
appearance the time-integrated pair annihilation spec- 
trum appears almost like the well-known 7-ray spectrum 
from decaying 7r°-mesons at rest, since it yields a broad 
bumpy feature located between 50 and 100 MeV. More- 
over, we demonstrated that for pair densities greater 10 10 
cm~ 3 in the jet the annihilation radiation will dominate 
the inverse Compton radiation, and due to its bump-like 
spectrum indeed may explain reported spectral bumps at 
MeV energies. The refined treatment of the inverse Comp- 
ton radiation including Klein-Nishina cross section cor- 
rections, the low-energy cutoff in the radiating relativis- 
tic particle distribution function and a realistic greybody 
accretion disk target photon spectrum leads to spectral 
breaks of the inverse Compton emission in the MeV en- 
ergy range with a change in spectral index Aa larger than 
0.5 as detected in PKS 0528+134 and 3C273. 
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